The electrostatic potential profile along a biased molecular wire: A model quantum 

mechanical calculation 
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We study the electrostatic potential of a molecular wire bridging two metallic electrodes in the 
limit of weak contacts. With the use of a tight-binding model including a fully three-dimensional 
treatment of the electrostatics of the molecular junction, the potential is shown to be poorly screened, 
dropping mostly along the entire molecule. In addition, we observe pronounced Friedel oscillations 
that can be related to the breaking of electron-hole symmetry. Our results are in semi-quantitative 
agreement with recent state-of-the-art ab initio calculations and point to the need of a three- 
dimensional treatment to properly capture the behavior of the electrostatic potential. Based on 
these results, current-voltage curves are calculated within the Landauer formalism. It is shown that 
Coulomb interaction partially compensates the localization of the charges induced by the electric 
field and consequently tends to suppress zones of negative differential resistance. 
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I. INTRODUCTION 

Due to important technical progress, the field of molec- 
ular electronics, born in the mid 70s with the proposal 
of Aviram and Ratner to use single organic molecules 
as rectifiers, El receives rapidly growing interest .Emu In- 
deed, new fabrication methods and probes now enable 
individual molecules or small numb,ejja_|Qf them to be 
connected to macroscopic electrodesBBQiHlj Among these 
methods, ^me may cite, for instance, the break-junction 
techniqueQ'Hj and the use of a conducting Atomic- Force- 
Microscope (AFM) to contact molecules absorbed on a 
metallic surface. 

On the theoretical side, the problem posed by these 
experimental works is highly challenging. We are facing 
a non-equilibrium many-body problem where, moreover, 
the coupling to a phonon bath may also be of impor- 
tance. Up to now, most of the studies have focussed 
on the coherent regime and the Landauer approach has 
been employed to obtain the conductance from ab initio 
or semi-empirical models.H Important inelastic. processes 
were included only within simple modelsoE30E3E3 and 
much further progress is needed before one may hope to 
reach a satisfactory understanding of the problem. 

The exact number of molecules contacted by the leads 
remains for a large P&ft uncontrolled in the experimental 
set-ups cited above.EHa In theoretical modelling it is con- 
venient to assume that a single organic molecule bridges 
two semi- infinite metallic electrodes (cf . Fig. . Another 
important experimental aspect is the fact that current- 
voltage characteristics are measured with applied volt- 
ages up to a few volts, values which bring us well away 
from the linear regime. 

In this context, a central question concerns the electro- 
static potential profile of a biased molecular wire. The 
importance of this-issue was first demonstrated by Datta 
and coworkers) 1 iri Using semi-empirical models, they 



have shown that different choices of the electrostatic pro- 
file have a profound effect on the current-voltage char- 
acteristics of a molecular junction. For instance, the 
transport properties are strongly modified depending on 
whether the potential drop occurs at the interface be- 
tween molecule and electrode or along the molecular wire. 
In fact, it is natural to assume that even the details of the 
potential shape have a considerable effect on molecular 
conductance. 

Recently,,— a^-,few works along that line have been 
r ep ort ed ,r 6 IF ll 1 9 They give us a rather ambiguous view 
of this fundamental problem. Model calculations involv- 
ing self-consistent solutions of the coupled Poisson and 
Schrodinger equations suggest that the potential drop 
occurs mainly a± the interface between the molecule and 
the electrodes.Eil Within the molecule, the electrostatic 
potential is then found to be essentially flat. Screen- 
ing appears to be very efficient within this approach and 
the final conclusions axe in agreement with some pre- 
vious investigations. oE£l However, these model calcula- 
tions involve a drastic approximation: instead of solving 
the full Poisson equation, the authors of Ref. 16 have used 
a one-dimensional version of it. Implicitly, it is then as- 
sumed that the lateral dimensions of the molecule are 
much larger than the screening length. For quasi-one 
dimensional systems with lateral dimensions of the or- 
der of a few angstroms, such as the organic molecules 
used in recent experimental work, this approximation is 
clearly questionable .La Indeed, recent state-of-the-art ab 
initio calculations-ntLcarbon and gold chains show a quite 
different picture. EZrll3 In these works, the potential drop 
occurs not only at the interface but rather along the en- 
tire molecule. Moreover, the local potential is found to 
display pronounced Friedel oscillations. Contrary to pre- 
vious results, screening appears to be rather inefficient, 
even for metallic wires. 

We are then left with two different pictures and it be- 
comes clear that a full understanding of the electrostatic 
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potential profile in biased molecular wires or metallic con- 
strictions is still lacking. In this work, we re-address the 
problem following the approach by Mujica et alS3 We 
perform model calculations and solve the coupled set of 
two equations: the Poisson equation for the electrostatics 
and the Schrodinger equation for the electronic structure. 
In .this respect, our work is similar to that of Mujica et 
al.i3 However, the calculation is modified in ways that 
we believe to be essential. In particular, we treat the 
real three-dimensional Poisson equation. As recently dis- 
cussed by usfS we expect that this proper handling of the 
electrostatic problem changes the qualitative behavior: 
The potential is poorly screened and falls off substan- 
tially along the molecule. Indeed, the electrostatic po- 
tential profile of the model calculations presented here is 
in semi-quantitative agreement with the ab initio results 
reported in Refs. [Tt] and [l^. However, while ab initio 
calculations are involved and intricate enough to leave 
the underlying physics essentially obscure, our model in- 
cludes only the ingredients necessary to capture the cor- 
rect screening effects, and consequently our calculations 
are rather economic in time. We thus believe that our ap- 
proach may help to gain further insight into the difficult 
problem of understanding transport through molecules, 
i.e. at a scale where quantum effects are prominent. A 
forthcoming work by Ghosh et alEi reaches conclusions 
similar to those presented here. 

In section |n|, our model Hamiltonian is introduced. 
The electronic density in the absence of a bias potential 
is then studied in section [II , using exact diagonalization 
and Hartree-Fock calculations that are shown to agree 
reasonably well which each other. The electrostatic pro- 



file is then calculated in section IV at the Hartree-Fock 
level. Finally, the current-voltage characteristics are dis- 
cussed in section [v|. 



II. MODEL HAMILTONIAN 

A. Coulomb interaction including image charges 

The physical problem posed by a molecular wire be- 
tween two infinite metallic reservoirs is far too compli- 
cated to be solved exactly and, to proceed, several ap- 
proximations are necessary. First, we assume that the 
surfaces of the two electrodes are infinite planes (cf. 
Fig. |l|). Second, the molecule is assumed to be weakly 
connected to the metallic electrodes so that their chemi- 
cal constitution is unimportant. This is certainly not ful- 
filled for some of the wires examined experimentally, with 
covalent bonds between molecule and electrode. Finally, 
we assume the characteristic time scale for electronic pro- 
cesses in the electrodes to be much shorter than the tran- 
sit time of electrons in the wire. The electrodes can then 
be treated as equipotential surfaces, and the Schrodinger 
equation is solved under these potential boundary con- 
ditions. We are then within the same framework used 
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FIG. 1: Idealized molecular junction. A molecular wire, mod- 
elled by a finite one-dimensional lattice, bridges two metallic 
electrodes with surfaces assumed to be infinite planes. The 
tunnelling contacts, effective only at the two molecular end 
sites indicated by 1 and N, are assumed to be weak. The 
chemical potentials in the left and right electrode are denoted 
by m and fi r , respectively. 



the Coulomb interaction potential which includes the im- 
age charges due to the metallic leads, keeping its three- 
dimensional character. Then, we solve the electronic wire 
problem. 

The Coulomb interaction energy reads 
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dVp(r)p(r>(r,r') 



(1) 



where p(r) is the charge density in the wire and ip(r, r') 
is the potential at point r = (x, y, z) produced by the 
charge located at point r' = (x',y',z'). The Coulomb 
potential is the solution of the Poisson equation 



A r v?(r,r') = -4nU8(r - r') 



(2) 



where for convenience we have introduced a factor U 
measuring the strength of the electron-electron interac- 
tion. In the absence of the metallic electrodes the solu- 
tion of (§) is given by the standard Coulomb potential 



yo(r.r') 
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(3) 



In the setup depicted in Fig. I], we require the potential 
in the absence of an external bias to vanish on the sur- 
faces of the metallic electrodes, i.e. at z = and z = L. 
A solution of the Poisson equation with these particular 
boundary conditions is found using the standard method 
of image charges 
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[(p (r + 2pLz, r') - tp (2pLz - r, r')] 



in Ref. 16 but proceed differently. We first determine 



(4) 

where L is the distance between the two electrodes (cf. 
Fig. [I]) and z the unit vector along the molecular axis. 

tp(r, r') is a genuine three-dimensional Coulomb poten- 
tial including the effects of the two semi-infinite metallic 
electrodes. The electrostatics of the molecular junction 
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is then governed by this potential. It remains to con- 
struct and solve the Schrodinger equation for the molec- 
ular wire. 



B. Tight-binding model including image charges 

In the following, we will mainly be concerned with 
conjugated molecules and, in particular, with their low- 
energy properties. An appropriate description can then 
be given by an effective tight-binding Hamiltonian for 
the 7r electrons only.Eilc3 In addition, since there exist ab 
initio results for short chains of gold atoms Jla it is also 
of interest to study systems with electrons in s-orbitals. 

We therefore attach Gaussian type orbitals of the form 



0„(r) = A s x s exp [ - a[x 2 + y 2 + (z - z n f]} (5) 

to each atomic site n. This allows us to model both 
s-orbitals with s — and p-orbitals with s = 1. The 
center of the orbital, z n = d + (n — l)a, depends on the 
distance d between each electrode and the molecule as 
well as the lattice constant a in the molecule (cf. Fig. [I]). 
The parameter a determines the spread of the state and 
gives an estimate of the electronic density. Finally, the 
normalization constants for s- and p-orbitals are given by 
A = (2a/n) 3 / 4 and A Y = 2{2/n) 3 / 4 a^ 4 , respectively. 

The explicit form (||) of the orbitals allows us to deter- 
mine the effective parameters entering the tight-binding 
model. In the following, we adopt the 'Zero Differential 
Overlap' (ZDO) approximationElE3E3 



0*(r)0 m (r) = \4> n (r)\ 2 5 n , T 



(6) 



which remains valid as long as the orbitals 4> n are strongly 
localized on the atomic sites n. It implies orthogonality 
between the Gaussian orbitals on different sites and, most 
importantly, leads to a drastic reduction of non- vanishing 
matrix elements of the Coulomb operator since only the 
two center integrals are retained in the final model. In 
particular, no exchange integrals will appear. These inte- 
grals, involving a differential overlap, are usually negligi- 
bly small compared to the Coulomb integrals. Moreover, 
we approximate the positive cores by point charges lo- 
calized at the atomic sites m with coordinates (0, 0, z m ). 
Then, the energy of an electron localized at site n due to 
all positive core charges and their images becomes 

N 

e n = - u J2 d 3 r\(f) n (r)\ 2 (p(r,z m z) . (7) 

m=l 

Within the ZDO approximation the only finite Coulomb 
matrix elements are related to the interaction energy be- 
tween electrons localized at sites n and n' 

U n , n , = U [d 3 r /dV|0„(r)|V(r,r')|0„'(r')| 2 . (8) 



It is worthwhile to notice that the interaction terms (0) 
and (H) depend on the position n along the chain due to 
the image charges but also because of the finite size of 
the molecular wire. 

Within the ZDO approximation kinetic energy contri- 
butions vanish as a consequence of (^|). To lowest order, 
the overlap leads to a constant shift of the on-site en- 
ergy e n that may be disregarded and to nearest neigh- 
bor hopping. The hopping matrix element t cannot be 
evaluated directly within the ZDO and, thejpfpse, has 
to be treated as a parameter of the model.EiKp^l How- 
ever, it is possible to relax the approximationEj and es- 
timate the hopping matrix elements from (||). Doing so, 
we have found that their dependences on n are not pro- 
nounced, and, moreover, this kind of more sophisticated 
treatment would not change qualitatively our final con- 
clusions. Therefore, we assume the hopping matrix el- 
ements, t, to be constant along the chain, and use the 
ratio U/t as a parameter to examine the importance of 
electron-electron interaction. 

With the parameters just discussed, we obtain a 
description of the electrons in terms of an effective 
tight-bindipg. model which includes long range Coulomb 
interactionEH 

H = ^2( e n + Vn)4l trr Cn,<T + t( c l+l,a C n,a + h.C.) 
n , a n , ct 

+ ^ ^ y ^n,n' C n,a C n' ,a' Cn ' ,o' c n,a ■ (9) 

c n a ( c n.a) & r e the usual creation (annihilation) opera- 
tors for an electron with spin a in the local state cj) n . 
In the first term, we have accounted for an additional 
shift of the local potential due to an external bias. With 
the chemical potentials fii and fi r in the left and right 
electrode, respectively, the shift at site n is given by 



Mi 



fj, r - \n 



(10) 



The resulting tight-binding model (g) is mostly de- 
fined by the geometry of the molecular junction. The 
only other parameters are the energies t and U, which 
determine the strength of the kinetic and Coulomb en- 
ergies, respectively. It is usually believed that con- 
jugated molecules lie in an intermediate regime where 

u/t i ... i.yy 



III. ELECTRON DENSITY WITHOUT BIAS 
POTENTIAL 

Consider first the situation without bias potential, 
lii = /i, r or v n = 0. We analyse for this case the elec- 
tron density at equilibrium in the absence of electron 
transfer between the molecule and the leads. In all the 
calculations presented here, we assume the chains to be 
electrically neutral with on average one electron per site. 
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FIG. 2: Electronic density without bias potential for a molec- 
ular wire of carbon (s = 1) with N = 12 sites. The param- 
eters are a = 4.5/a 2 and U = t. The full curve is the exact 
result for a molecule-electrode distance of d = 2a. The dotted 
curve is for the same geometry but at Hartree-Fock level. The 
dashed curve is the exact result for the case without metallic 
electrodes, d — ► oo. 

Because of the electron spin, this corresponds to the half- 
filled case. Furthermore, we restrict ourselves to the 
ground state of the system. 

We have done exact diagonalization studies for chains 
of up to 12 sites. Results for the charge density of a 
chain with 12 sites shown in Fig. are typical for other 
cases studied. Two main features can easily be seen. 

(i) The electron density is non-uniform: the electrons 
have a tendency to shift towards the middle of the chain. 

(ii) Because of this non-uniformity, substantial Friedel 
oscillations occur across the wire. These features can 
be explained by invoking electron-hole symmetry as it is 
shown below. 

For the case of a half-filled band it has long been known 
that the electron density is uniform, i.e. does not depend 
on the site index n, for models defined on bipartite lat- 
tices in such a way that electron-hole symmetry is ful- 
filled. This theorem was first discovered-fpr the ground 
state of the free-electron (Hiickel) mode£-3 and later ex- 
tended to some interacting systems.o More recently, a 
generalization to canonical and grand canonical ensem- 
bles was proven for a large class of modelsEj Mathe- 
matically, the theorem applies to models invariant with 
respect to the transformation <y n a — > (— l)"c„ i(T . 

Applying this transformation, we find that the Hamil- 
tonian of the molecular wire (|^) is invariant only if the 
equality 

e n = -\j^U n , m (l- 5 -^)+K (11) 

m=l ^ ' 

is fulfilled where K is a, constant. This implies that 
the electron-ion interaction is essentially compensated, 
up to a constant term, by the repulsive electron-electron 
interaction. Notice that allowing for hopping matrix 
elements that depend on the position in the chain 



would not change this equality. Indeed, a term like 

J2r h a *n,n+i(4+i,<7 C «,<r + 4,<r c n+i,o-) remains unchanged 
when the electron-hole transformation is applied. From 
this point of view, it is not necessary to go beyond the 
ZDO approximation. 

The equality ([ll]) can be satisfied only in very partic- 
ular cases where at least one of the following conditions 
is satisfied: 

• the chain is infinite; 

• there is no interaction, U = 0; 

• a — > oo, corresponding to interactions between 
point charges. 

None of these criteria are fulfilled in realistic cases of 
interest. With the exception of carbon nanotubes, ex- 
perimental molecular junctions involve relatively short 
(< 10 nm) molecular chains.0 The Coulomb interaatioH. 
is of the same order of magnitude as the kinetic termseMij 
and certainly not negligible. Even then, the electron den- 
sity would still be uniform if the electron-orbitals are re- 
duced to points (a — > oo) irrespective of the presence 
of image charges. Therefore, one may say that the non- 
uniformity of the density comes from the lateral extension 
of the electronic clouds, the p- or s-states. This shows, 
once again, the need of a three-dimensional treatment of 
the electronic structure. 

The exact shape of the electron density is determined 
by (i) the strength of the Coulomb interaction, U/t, (ii) 
the spread of the electronic orbitals, a, and (iii) the geom- 
etry of the system, expressed in our model by the length 
(n — l)a of the molecule and by the distance d between 
the electrodes and the molecule. 

All these factors contribute to yield, instead of the uni- 
form density condition ([□]), the relation 

^n = -\jZUn, m U- 5 ^f) + fdn) (12) 
m=l ^ ' 

where is a function of the position in the chain; its 

dependence on n is responsible for the non-homogeneity 
of the electronic density. It is of interest to understand 
the effects of each of these parameters separately. 

Coulomb interaction. Starting from the non- 
interacting case where the electronic density is uniform, 
increasing U results in an increase of the Friedel oscil- 
lations until the electrons start to localize. In the limit 
U — > oo, half filling leads to a Mott insulator and one 
recovers a constant electron density. 

Spread of the electronic orbitals. When the orbitals 
are reduced to points (point charge limit), a — > oo, the 
electronic density is uniform. Indeed, in this particular 
limit the electron-ion and electron-electron interactions 
have the same form restoring the electron-hole symmetry 
of the molecular Hamiltonian. Accounting for a spread of 
the orbitals increases slowly the amplitude of the Friedel 
oscillations. 
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Geometry of the molecular junction. Two effects need 
to be distinguished: (i) the finite size effects and (ii) the 
image-charge effects. In Fig. |[ the dashed curve shows 
the density for the very same system used to obtain the 
other two curves except that now no metallic electrodes 
are present. Therefore, the dashed curve contains only 
the finite size effects. From this particular example, one 
sees that both effects contribute to the Friedel oscillations 
and that, in order to get the correct density, image inter- 
actions should be included unless the electrode-molecule 
distances is large. In fact, in the example of Fig. || the 
contribution to the oscillations due to the image charges 
is the larger one. The Coulomb interaction with the elec- 
trodes described by the image charges is therefore crucial 
to estimate transport properties of molecular wires or 
metallic constrictions. We have observed, as expected, 
that finite size effects alone tend to disappear when the 
chain size is increased. In contrast, the effect of image 
charges tends to become more important: the amplitude 
of the Friedel oscillations increases with the system size 
due to the presence of the two metallic electrodes. This 
tendency should continue until for long wires, that we do 
not consider here, the oscillations occur predominantly 
near the edges. 

Finally, the dotted curve of Fig. shows the electron 
density with metallic electrodes calculated within the 
Hartree-Fock approximation; it is in very good agreement 
with the exact result. We have performed calculations for 
different values of U up to U = it and observed that, as 
far as the electron density is concerned, the Hartree-Fock 
approximation gives reasonable results. This allow us to 
use this mean-field approach for the study of the electro- 
static profile in biased molecular wires (section 
their transport properties (section^). 



IV. ELECTROSTATIC POTENTIAL IN BIASED 
MOLECULAR WIRES 

The application of the screened on-site electrostatic po- 
tential to calculations of conductance properties is useful 
only within a single electron theory, e.g. in the mean-field 
description of the process. This is how this_concept was 
applied in the calculatio n o£ Mujica et alx3 and in the 
density functional theory .fLZMiS Here, our calculations are 
done using the Hartree-Fock approximation applied to 
the molecular Hamiltonian (^|). In the previous section, 
by comparing the mean-field electron density with the 
exact one, we have already shown, that the approximate 
result is reasonable in the range of parameters proper to 
conjugated molecules. More precisely, the Hartree-Fock 
density is in very good agreement with the exact results 
for values of U up to, approximately, It (cf. Fig. ||) and 
remains of reasonable accuracy up to values of about At. 
In the following we present only results for U — t, but 
the same qualitative picture arises also for larger values 
off/. 

Starting from the initial tight-binding model (^), we 



build an effective Hamiltonian by salving self-consistently 
the usual Hartree-Fock equationsEJ leading to 



IV) and 



t n,n'{V)(cl ta C nitT +h.C.) 



(13) 



where V = fii — \i r . The chemical potentials will al- 
ways be chosen such that \xi — — fi r . c n (V) is the effec- 
tive on-site potential which includes the ionic attraction 
and the electron-electron repulsion incorporated within 
a mean field picture as well as the on-site potential ([h]) . 
t n ,n'(V) is the effective hopping matrix element which in- 
cludes the exchange terms. Note that a vanishing of ex- 
change integrals in the ZDO approximation is not equiv- 
alent to an Hartree approximation. Indeed, a mean field 
approximation of (^|) still contains exchange terms of 
the form - jE^n'VfiA',^!',^,' which, due 
to the long range part of the Coulomb potential (Q), in- 
clude long range hopping. 

Without Coulomb interactions, the electrostatic po- 
tential is given by the ramp defined in (|l^), i.e. by the 
potential in the absence of a molecule. This linear pro- 
file has been used sometimes in the liter^uxe to study 
non-linear current- voltage characteristics. lZtl3 

In the presence of Coulomb interaction, the linear pro- 
file is changed and the screened electrostatic potential, 
E n (V), is given by the difference between the on-site term 
with and without bias voltage. 



E n (V) = e n {V) - e„(0) 



(14) 



A typical example is shown in Fig. |3| for a chain with 
20 sites. We see two main features: (i) screening is not 
very efficient and, consequently, there is only a small po- 
tential drop at the interfaces but a finite slope of the 
potential along the entire molecule; (ii) there are sub- 
stantial Friedel oscillations along the profile. _ 

Our results differ from those of Mujica et alS3 despite 
the fact that the two models are closely related. These 
differences mostly come from the fact that we solve the 
Poisson equation without resorting to a one-dimensional 
approximation. As recently shown by us within a classi- 
cal model,E3 a three-dimensional treatment of the electro- 
statics is necessary, and in fact leads to the identification 
of the lateral thickness of the molecule as a new generic 
attribute that determines the potential profile. In Fig. |^, 
we show a comparison between our results and the cal- 
culations of Damle et alxB for a chain of six gold atoms. 
It is important to stress that we did not try to fit the ab 
initio curve but, instead, we simply chose a reasonable 
set of parameters. Note, however, that we include the 
response of infinite leads while the ab initio calculations 
take only small metal clusters into account. Furthermore, 
the asymmetry in the latter case indicates that the ab ini- 
tio calculations result in a charged molecule while in our 
case the molecular wire always remains neutral. Other 
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FIG. 3: Electrostatic potential profile for a carbon wire (s = I 
in Eq. |^) with N = 20 sites. The parameters are a = 4.5/a 2 , 
U = t, V = Hi — fj, r = t and d = 2a. The dashed curve is 
the unscreened potential without molecule (ramp potential). 
The full curve is the screened potential in the presence of the 
molecular wire. It shows a small decrease in slope along the 
entire molecule with substantial Friedel oscillations. 



ab initio calculationsEIIE£rE2l on similar models are also in 
agreement with our observations. 

To conclude this section, we summarize our main re- 
sults: in the relevant range of parameters, screening is not 
very efficient in the wire and the drop of the potential oc- 
curs along the entire molecule. Additionally, substantial 
Friedel oscillations are present in the electrostatic profile. 
Our respite are in good agreement with recent ab initio 
results. LL3li3 When compared with the results of Ref. [if], 
our treatment stresses the need to study the full three- 
dimensional problem, as far as the electrostatics of the 
molecular junction is concerned. 

In the next section, we discuss consequences of the 
screening effects for the current- voltage characteristics of 
molecular wires studied within the Landauer formalism. 



V. CURRENT- VOLTAGE CHARACTERISTICS 
OF A MOLECULAR WIRE IN THE WEAK 
TUNNELLING CONTACT LIMIT 

In a model where only the first and the last atom of a 
molecular chain couple to the corresponding metal leads 
and at T — 0, the Landauer conduction formula yieldtoEJ 



2e 



alE\G 1N (E,V)\ 2 A l (E,V)A r (E,V) (15) 



where V = \i\ — \i r and A;/ r are the spectral functions 
for the left and right reservoirs. The molecule-lead cou- 
pling only occurs at sites 1 and N and Gin is the ma- 
trix clement of the exact Green function of the molecular 
junction between these sites (cf. Fig. |l|). 



FIG. 4: The electrostatic potential profile for a gold wire 
(s = in Eq. |) with N = 6 sites obtained from Hartree-Fock 
calculations is shown by the filled circles. The parameters 
used are d = 0.9a, a — 4.5/a 2 , and U = t. For comparison 
the full curve shows the ab initio results taken from Ref. 



tacts, where it is possible to use second order perturba- 
tion theory in the tunnelling matrix element .EH In the 
limit of 'extremely' bad contacts of interest here, the 
Green function of the system in formula ( |l5| ) may be 
replaced by tbe. Green function of the isolated molecu- 
lar wire, Gi N S$ Moreover, assuming that the product of 
spectral densities does not significantly depend on energy 
in the range between fii and /i r , one gets 



T 2 e a2 
irn 



dE\C% N (E,V)\< 



(16) 



This equation can 
the Landauer forrnuh 



understood as a special case of 
adapted to the case of bad con- 



where Ao is the spectral density at zero bias. 

From this simplified equation, we can calculate the 
current- voltage characteristics of the molecular junction, 
using only the Green function of the isolated molecular 
wire, evaluated in the presence of the electrostatic po- 
tential created by the metallic electrodes. Examples are 
given m Fig. | for increasing Coulomb strength U starting 
from the non-interacting case where the on-site potential 
is given by the ramp (|l0|). 

In all cases, the I-V curves have a staircase struc- 
ture -ffibirJa.. is a common feature in the weak tunnelling 
limitli3'LjO and simply reflects the discreteness of the 
molecular electronic spectrum. Indeed, an increase of 
the bias potential corresponds to an increase of the win- 
dow of integration in formula (|l6|). Therefore, a jump in 
the I-V curves means that one more discrete molecular 
level enters this window of integration. 

It is interesting to note that Fig. || shows also wide 
regions of negative differential resistance, in particular in 
the non-interacting case. Within our simple formulation, 
they can be explained by the localization of charges in- 
duced by a strong enough electric field: at sufficiently 
high bias voltage, charge carriers are localized at one end 
of the chain resulting in a decrease of the current. This 
could be an artifact of our weak molecule-electrode cou- 
pling model, however negative differential resistance has 
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FIG. 5: Current-voltage characteristics for a carbon molecule 
(s = 1) with N = 12, d = 2a, a = 4.5/a 2 , and for different 
Coulomb interaction strengths C//t = 0, 0.1, 1 are depicted by 
the dotted, dashed, and full line, respectively. The current is 
given in units of (2e/7r7i)A(j/f. 



been found in recent experiments £J and electric field in- 
duced localization could be a way to understand these 
experimental findings. 

The Coulomb interaction has two main effects on the 
I-V characteristics. First, the positions of the current- 
voltage steps are shifted to higher voltages reflecting the 
displacement of the molecular levels to higher energy 
with increasing Coulomb interaction. Second, the lo- 
calization of the charges due to a strong electric field is 
partially compensated by the electron-electron repulsion. 
These screening effects attenuate the negative differential 
resistance effects (cf. Fig. . These observations are sim- 
ilar to those made by Mujica et al.,KZl where a Hubbard 
model was studied at the Hartree level. 



VI. CONCLUSION 

We have addressed the problem of calculating trans- 
port properties of a molecular wire bridging two semi- 
infinite metallic electrodes. A first important part of 
this task is to determine the^ab^tajstatic potential profile 
through the biased wireoSE^EZrlia Indeed, it is of im- 
portance to know how screening effects modify the ramp 
potential (Eq. 10) existing in the absence of the organic 
molecule. _ 

This work and our earlier papeiO resolve discrepan- 
cies between answers available in the literature. On 
the one hand, a tight-binding model combined with a 
one-dimensional Poisson equation gives a strong screen- 
ing version of the problem: the drop of the potential 
occurs at the interfaces and the— potential is therefore 
almost flat within the molccule.E£l On the other hand, 
recent ab initio results give a weak screening version: 
no substantial drop at the interfaces but rather a de- 
crease along the entire molecule together with substantial 



Friedel oscillations 0-0 

In this paper, we have proposed a modified tight- 
binding model to address this question. It is based 
on three main ingredients, (i) We introduce a three- 
dimensional Coulomb potential (Eq. ||) which includes 
the image interaction with the two metallic electrodes 
assumed to have planar surfaces, (ii) The electrons lo- 
calized at atomic sites are modelled by Gaussian type 
orbitals of finite lateral extent (Eq. ||). (iii) The positive 
background is assumed to be a set of point charges local- 
ized on the atomic sites. With these three ingredients, 
it is possible to evaluate the various terms of our model 
Hamiltonian: the on-site energy (Eq. and the electron- 
electron interaction (Eq. ||). They are all functions of the 
position on the chain due to the finite size of the system 
and the image charges induced by the electrodes. All 
our calculations are done in the weak tunnelling limit, 
assuming bad contacts. 

This model yields an electronic density that is non- 
uniform already in the absence of a bias potential — the 
electrons prefer to be in the center of the wire — and 
displays pronounced Friedel oscillations (Fig. ||). These 
characteristics are due to the fact that our modelnjdoea 
not fulfill, in general, the electron-hole symmetry .E-3E§E£l 
In the presence of an applied voltage the electrostatic 
potential profile does not drop appreciably at the inter- 
faces but rather, in accordance with the ab initio results 
of Refs. [It] and |l8|, it decreases along the molecule with 
substantial Friedel oscillations appearing along the entire 
profile (Figs. p| a nd [|). Our results are different from the 
ones of Ref. Hq where a one-dimensional Poisson equa- 
tion was used. This disagreement stresses the need to 
perform a three-dimensional calculation, as done here, 
to properly describe the electrostatic properties of the 
molecular junction. Finally, the current-voltage charac- 
teristics shows a staircase struct-use (Fig. |J) as is common 
in the weak tunnelling limit .ESEZI Zones of negative dif- 
ferential resistance are found due to charge localization 
induced by the electric field. The main effects of the 
Coulomb interaction, within our approximations, are, on 
the one hand, a shift to higher energies of the position 
of the current-voltage steps and, on the other hand, a 
partial compensation of the localization of the electrons 
diminishing the negative differential resistance effects in 
agreement with a previous studyB3 

In closing, it is important to stress some limitations of 
our model. On the one hand, we consider coherent trans- 
port, assuming that the electrons are transferred from 
one lead to the other in a single quantum mechanical 
process. This is a good approximation if the tunnelling 
time is much less than the inelastic scattering time. For 
organic molecules, this transit time could be of the same 
order of magnitude as the intramolecular vibronic re- 
laxation time, especially in the weak contact limit em- 
ployed here.U In this case, part of the current could be 
due to sequential tunnelling, where the molecular wire 
would be successively charged and discharged. This im- 
portant issue remains to be studied further. We have 
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neglected charging effects assuming the molecule to re- 
main-neutral. At high voltage, this approximation could 
fail£3 The average charge number of the molecule could 
increase in analogy with Coulomb blockade phenomena 
observed in mesoscopic metallic double tunnel junctions 
and quantum dotsjp and, more recently, in multiwall 
carbon nanotubesHU 

A proper handling of the full problem requires to treat 
the wire as an open system dynamically coupled to the 
electrodes and to the vibronic degrees of freedom. This 
program is far beyond the scope of the present work. 
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